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We study the classical ground states of the exchange-coupled Heisenberg antiferromagnet on the Pyrochlore 
lattice, a non-Bravais lattice made of corner-sharing tetrahedra. In particular, we map out the entire phase 
diagram for the case of first and second nearest neighbor interactions. In this phase diagram we find four 
complex non-coplanar ground states based on different ordering modes. These are the Cuboctahedral Stack 
state, a (111) stacking of Cuboctahedral states, and three families of commensurate spiral: the Kawamura 
states, constructed from three different combinations of {§ §0} modes, the Double-Twist state, also constructed 
from {§§0} modes, and the Multiply-Modulated Commensurate Spiral state, constructed from { § | §} modes. 
We also briefly look at states involving the two kinds of third nearest neighbor interactions on the Pyrochlore 
lattice. In this region of parameter space we again find the Cuboctahedral Stack state, and we also find another 
non-coplanar state in the form of a new kind of Alternating Conic Spiral. 

PACS numbers: 75.25.-j, 75.30.Kz, 75.10.Hk, 75.40.Mg 



I. INTRODUCTION 

This paper concerns the classical ground states of the 
Heisenberg Hamiltonian on the Pyrochlore lattice, a non- 
Bravais lattice made of corner-sharing tetrahedra. In this pa- 
per we focus our analysis on the case of first and second near- 
est neighbor interactions J\ and J 2 , and we briefly mention 
some results that include the effects of third nearest neigh- 
bor interactions J3 and J3. (Each site on the Pyrochlore lat- 
tice has two distinct kinds of third nearest neighbor, at the 
same separation, yet inequivalent by symmetry - having dis- 
tinct neighbors at the same distances appears to be universal 
in non-Bravais lattices.) 

An important motivation to work out the phase diagram of 
magnetic systems is the inverse problem arising from neu- 
tron diffraction: not the crystallographic problem of fitting the 
magnetic structure from the measured Bragg intensities, but 
rather the physical problem of inferring - in the absence of 
magnon dispersion relations that would require dynamic neu- 
tron scattering - for what possible interactions is the observed 
structure the ground state. (One motivation for this project 
was the puzzle presented in Ref. iflOl in GeNi2C>4, where it 
was necessary to posit interactions as far as J4.) 

The other motivation is to facilitate the theorists' pursuit of 
model systems with unusual properties. For example, when- 
ever frustrated magnets have a large number of degenerate 
ground states, quantum or thermal fluctuations, or dilution, 
can select a particular ground state - a form of so-called "or- 
dering due to disorder" 10). Thus, one interest in our study is 
to identify parameter combinations on the phase diagram that 
lead to such degeneracies. This tends to happen along "degen- 
erate phase boundaries" (see Ref. Q] sections II C and VII): 
their property is that one gets two limiting states depending 
on which side of the boundary the limit is taken from, yet it is 
possible to continuously turn one of these into the other within 
the manifold of degenerate states that exists on the boundary. 

Another special property, which we are particularly con- 
cerned with in this paper, is inherently non-coplanar states. 
These can produce an anomalous Hall effect due to Berry 



phases of itinerant electrons lfTT14T4ll : they open up possibili- 
ties of multiferroic behavior ll5l[T8l[T9lk they turn into chiral 
spin liquids if spin order is destroyed; and their symmetry- 
breaking leads to an order parameter in the form of an 0(3) 
matrix, which allows unusual topological defects (Z2 vortex 
and Z2 domain wall |5 1). 

The Pyrochlore spin lattice represents many antiferromag- 
nets in any of three well-known crystal structures: the B sites 
of the spinel structure, one cation sublattice (of the two equiv- 
alent ones) in the Pyrochlore structure; or the sublattice of 
small atoms in the (metallic) Laves phase structure. It is 
highly degenerate when only the nearest-neighbor interaction 
is included; in studies where other terms were added that 
break the degeneracy, surprisingly few of them added farther- 
neighbor coupling (apart from dipolar spin ice 11281 ). Instead, 
anisotropies were added such as local easy planes ||29ll or 
Dzyaloshinskii-Moriya anisotropic exchange l30l . Alterna- 
tively, lattice distortions were considered BTH3411 . 

In the case of the Pyrochlore, Ref. [24] worked out (in ef- 
fect) the optimum Luttinger-Tisza modes for general Hamil- 
tonians with Ji up to J4, but limited to J 3 = J3. The 
main systematic exploration to date has been of the ( Ji, J2) 
Pyrochlore [or the (Ji, J3) case, which is equivalent to the 
{Ji, J2) case so long as J\ is large and antiferromagnetic in 
sign]. It uncovered a multiple-wavevector, not quite commen- 
surate state with especially soft fluctuations [8,9 36 , 37) that 
is not quite fully understood. Another context in which farther 
neighbor couplings naturally arise is in metallic Pyrochlores, 
in the form of electron-mediated R.K.K.Y. exchange inter- 
actions (that oscillate and decay with separation as a power 
law). 



II. NOTATIONS AND METHODS 

In this section we describe the Pyrochlore lattice and we 
also explain the methods and diagnostics we used to gener- 
ate and analyze the spin configurations. We present a num- 
ber of ways of visualizing the Pyrochlore lattice, as certain 
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ways of thinking about the structure of the Pyrochlore lattice 
are helpful for understanding certain spin configurations. We 
also discuss the Iterative Minimization algorithm, our main 
method of generating spin configurations. Finally we discuss 
a number of ways of understanding and interpreting the re- 
sults of the Iterative Minimization simulation, including the 
Luttinger-Tisza method, numerical Fourier transforms, least 
square fits and Variational optimization. 

A. Pyrochlore Lattice 

The Pyrochlore lattice consists of four face-centered cubic 
(FCC.) sublattices, numbered 0, 1, 2 and 3. We take the edge- 
length of the cubic unit cell of each FCC. sublattice to be 
one. We choose the origins of these four FCC sublattices 
to be located at (0, 0, 0), (0, |, ±), (±,0, \) and (A 1,0), re- 
spectively. 

Alternately, the Pyrochlore lattice can be viewed as an ABC 
stacking of triangular and Kagome lattice layers along a (111) 
direction in real space as in figure [T] 

Our Hamiltonian is 

1 N 

H = --22z2 J H s i- s 3 i (2- 1 ) 

The s, are classical 3-component unit vectors which sit on the 
sites Yi of the Pyrochlore lattice. The lattice contains N sites. 
A positive (negative) Jy represents a ferromagnetic (antifer- 
romagnetic) interaction between the spins Sj and Sj. Since 
the interaction Jjy is a function of the displacement vector 
Rjj — Yi Yj between the spins s, and Sj, we can write 

Jij = J[Rij). 

Each site on the Pyrochlore lattice has six first nearest 
neighbors, 12 second nearest neighbors and 12 third nearest 



FIG. 1 : The Pyrochlore lattice can be viewed as a stacking of two- 
dimensional triangular and Kagome lattice layers along a (111) axis. 
The picture shows one (black) Kagome lattice layer sandwiched be- 
tween an upper (red) triangular lattice layer and a lower (blue) trian- 
gular lattice layer. 




neighbors, but the third nearest neighbors come in two kinds, 
and each site has six of each kind. First nearest neighbors 
have an interaction Ji, second nearest neighbors have an in- 
teraction J r 2, third nearest neighbors of the first kind have an 
interaction J3 and third nearest neighbors of the second kind 
have an interaction J3. From now on we refer to first nearest 
neighbors as "Ji neighbors", second nearest neighbors as " J 2 
neighbors", etc. 

The J \ neighbors of a site can be reached by traveling along 
one edge of a tetrahedron. The J2 neighbors of a site can be 
reached by traveling along two non-collinear edges of two ad- 
jacent tetrahedra. The J3 neighbors of a site can be reached by 
traveling along two collinear edges of two adjacent tetrahedra. 
Finally, the J' z neighbors of a site can be reached by hopping 
across a hexagon in the Kagome lattice layers. The two kinds 
of third nearest neighbors are obviously inequivalent by sym- 
metry since every two J3 neighbors have a spin lying halfway 
between them on the line connecting them, while there is no 
spin lying on the line connecting two J3 neighbors. 

In this rest of this paper we limit our study to J\ and J2, 
except that J3 and J3 are used in Section [V] and Appendix [b] 
and J 3 is mentioned in Section [TV C| 

B. Luttinger-Tisza Analysis 

The Luttinger-Tisza method [3 6| is a method for finding 
the ordering wavevector that characterizes the lowest energy 
state of ( |2.1| i for a given set of interactions J^. Recall that the 
goal of any method (analytical or computational) for finding 
the ground states of ( |2.1[ ) is to find the set of N spins {s^} 
which minimize the total energy subject to the constraint that 
all of the spins have unit length, i.e. |sj| 2 = 1 V i. The con- 
straint that all of the spins have unit length is often called the 
strong constraint. 

The Luttinger-Tisza method attempts to accomplish this by 
instead trying to minimize ( |2.1| i subject to what is known as 
the weak constraint, 

Y,h\ 2 = N. (2.2) 

i 

By formulating the problem in terms of the Fourier Trans- 
forms of the spin configuration in each Bravais sublattice, one 
can obtain a lower bound on the energy of the true ground 
state. One then searches for a solution to the weak constraint 
problem which is composed solely of modes related by the 
symmetry of the underlying Bravais lattice. In an ideal situa- 
tion we would be able to find a solution to the weak constraint 
problem that also satisfies the strong constraint, is composed 
solely of modes related by the symmetry of the underlying 
Bravais lattice, and has an energy equal to the lower bound. 
If such a configuration could be found, it would rigorously be 
the ground state. 

To formulate the method we first rewrite equation \2.\\ in 
a way which clearly shows which spins are in which Bravais 
sublattice. We use the notation a(i) to indicate which Bravais 
sublattice spin i belongs to, so the possible values for a(i) 
are 0, 1, 2 and 3 in our case. When it is convenient we will 
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also write i £ a to indicate that spin Sj is located in the a 
Bravais sublattice. From here on we specialize to the case of 
the Pyrochlore lattice, which has four EC. C. sublattices. Then 
we can write 

Jij = Ja(i)P(j)(R<ij) (2.3) 

where J Q ^(R) is the interaction between spins in F.C.C. sub- 
lattices a and /3 that are separated by a vector R. 

Next we define the Fourier Transform of the spin configu- 
ration in the a th Bravais sublattice to be 

MqH^E^s (2-4) 

where q is a wavevector in the Brillouin zone of the F.C.C. lat- 
tice and N/4 is the number of spins in each F.C.C. sublattice. 
The inverse transform is 

s ' = ^lE eiq ' ri s«w(q)- ( 2 -5) 

Substituting this into ( |2.3[ ) yields the result 

H = -^^J^(q)S a (q).S /J (-q). (2.6) 

q a,P 

The coefficients of the quadratic form ( |2.6| l form a 4 x 4 
matrix J(q) with matrix elements 

Let A 1/ (q) be the four eigenvalues of this matrix at wavevector 
q, each with a normalized eigenvector u v a (q), where a runs 
over the four sublattices. These modes, transformed back to 
the actual lattice as u" = u v (ri), are exactly analogous to the 
Bloch wavefunction for the v band. Express the spin configu- 
ration in terms of these eigenmodes: 

S a (q) = £;w„i£(q), (2.8) 

so the coefficients have three Carteisan components for spin. 
The total energy is then 

H = -]T5> y (q)|w«,(q)| 2 . (2.9) 
q v 

Let X max (q) be largest eigenvalue of the matrix J(q) and 
let Ql.t. be a wavevector that maximizes X max (q) (gener- 
ically there is a star of symmetry-related wavevectors with 
magnitude \Ql.t. | that all do this). The quantity Al.t. = 
XmaxiQL.T.) is called the optimal Luttinger-Tisza eigen- 
value, and Ql.t. is the optimal Luttinger-Tisza wavevector. 
We refer to J(Ql.t. ) as the Luttinger-Tisza matrix. 

In Fourier space the weak constraint ( |2.2| i becomes 
S q v |wi/(q)| 2 = N so the ground state energy manifestly 
satisfies the inequality 

H > -NXl.t. ■ (2.10) 



A necessary condition for achieving this minimum is that 
the three components of S a (q) be linear combinations of the 
Luttinger-Tisza eigenmodes w^(Ql.t.)- The Luttinger-Tisza 
method works if one can find a normalized spin configura- 
tion composed solely of eigenmodes of Ql.t. (and symmetry- 
related wavevectors) which also achieves the lower bound on 
the energy, —NXl.t. - On non-Bravais lattices, it is not always 
possible to construct a normalized state using only Luttinger- 
Tisza eigenmodes. Many states require admixtures of a set of 
suboptimal modes in order to achieve normalization (TJ. Our 
experience working with the Pyrochlore lattice confirms that 
this is true. 

The key difficulty with the Luttinger-Tisza method is to en- 
sure | Si | = 1 for every site, when that is not generally the 
case for a Luttinger-Tisza eigenmode (expressed in real space, 
where it appears as a plane wave with a certain weight on each 
sublattice). On Bravais lattices, that is the case if the eigen- 
modes are written as complex plane waves exp(iQL.T. ■ r); 
from this we can always construct a ground state in the form 
of a coplanar spiral by using one component for the sine and 
another for the cosine |3|. Since the optimal Luttinger-Tisza 
wavevectors typically occur in symmetry-related stars, a non- 
trivial degeneracy would be possible by taking linear com- 
binations of two modes with different Q's from the same 
star; however, since each mode requires two spin components, 
that is not possible with physical three-component spins. But 
whenever Ql.t. is half of a reciprocal lattice vector (and thus 
lies on a special point of the Brillouin zone boundary), we 
have exp(jQL.T. • r) = ±1. Thus the plane-wave modes are 
real in this case, and we can combine up to three of them, 
associated with orthogonal spin directions, to construct con- 
tinuous families of degenerate classical ground states. 

For each set of couplings Jj,-, we can calculate the op- 
timal Luttinger-Tisza wavevector Ql.t.- This allows us to 
create a "Luttinger-Tisza phase diagram", showing in what 
domains of parameter space each kind of L.T mode is opti- 
mal. The L.T. phase diagram tells us the expected wavevec- 
tor content of the ground state for a given set of couplings, 
but cannot tell how those modes are combined to construct a 
normalized state, whether there are multiple degenerate ways 
to do so, or whether the actual ground state requires admix- 
tures of suboptimal modes for normalization. So although the 
Luttinger-Tisza phase boundaries typically correspond closely 
to the real phase boundaries, the true phase diagram often 
has extra subdivisions that the LT phase diagram lacks, e.g. 
the three Kawamura states (which all use the same Luttinger- 
Tisza modes, with admixtures of suboptimal modes). 



C. Iterative Minimization method 

Following Ref. 1 our main method for finding the ground 
states of ( |2.1| i is the Iterative Minimization simulation, a nu- 
merical approach which starts with a random spin configu- 
ration and generates states with progressively lower energy, 
until the method converges to some (local) energy minimum. 
Let us first define the local field to be minus the gradient of 
\2.\\ with respect to the components of spin i (analog of the 



4 



force in a mechanical system): 
SH 

h, = — 



(2.11) 



In any ground state, every spin must be aligned parallel to its 
local field. (If any spin were not, reorienting it along its local 
field would immediately lower the energy.) 

The method works in the following way. We start with N 
spins on the Pyrochlore lattice, all pointing in random direc- 
tions. Next we run a large loop, and on each iteration of the 
loop we pick N spins at random (one at a time) and reorient 
them so that they point along their local field. 

In our simulations we mostly work on a lattice that is an 
L x x L y x L z block of Pyrochlore lattice unit cells (where L x , 
L y and L z are integers) with periodic boundary conditions. 
We also have the capability to use twist boundary conditions, 
in which spins near a certain boundary of the lattice see their 
neighbors beyond that boundary twisted by an angle 9 about a 
certain axis in spin space. 

Once we have obtained a spin configuration using Iterative 
Minimization, we try to get a sense of the basic structure of 
the state using a few simple diagnostics. These are the Site 
Energy diagnostic, which gives a rough sense of the layout of 
the spin configuration in real space, the Common-Origin plot, 
which indicates which directions the spins are pointing in in 
spin space, and numerical Fourier Transforming, which gives 
the wavevector content of the state. Finally, we use Least 
Square fits of the spin configuration to the wavevectors ap- 
pearing in the numerical Fourier Transform to get an approxi- 
mate parameterization for the ground state. For simpler states, 
we can then use the technique of Variational Optimization to 
determine the idealized values of the parameters appearing in 
our approximate parameterization of the state. (All of these 
except the site energy were introduced in Ref.Q]) 



I. Site energy diagnostic 
We define the site energy of the spin S j to be 

Si = -~y] JijSt ■ Sj . (2.12) 

The average of this quantity over all the spins in the lattice is 
the "energy per site" and we denote it by £ = i S j. The 
total energy of the system can then be written as H = N£. 

We often find that spin configurations break up into sublat- 
tices in real space according to the site energy of the spins in 
those sublattices. Sometimes these sublattices are the famil- 
iar four FCC. sublattices of the Pyrochlore lattice, as in the 
Cuboctahedral Stack state (see section [V]), but sometimes the 
spin configuration breaks up into a more exotic set of sublat- 



2. Fourier Transforming 

Once we have found a spin configuration using Iterative 
Minimization, we take Fourier Transforms of that spin config- 
uration in each FCC. sublattice. Next we compute the Spin 
Structure Factor in each FCC sublattice to determine which 
wavevectors the configuration mostly consists of. The Fourier 
Transform of the spin configuration in the a th FCC sublat- 
tice is 



S«(q) 



1 



VP 



E 



s.,e 



(2.13) 



where q is a wavevector in the Brillouin zone of the FCC 
lattice and the sum is taken over all spins Sj in the a th FCC 
sublattice. In practice we compute these Fourier Transforms 
one spin component at a time. 
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We then examine 



S a q) 



tices, as in the Kawamura Sextuplet-q state (see section VI A i 



for each FCC sublattice (here 

the square includes a sum over the three Cartesian compo- 
nents from the spins, as well as real/imaginary parts from 
the Fourier transform.) Sorting these in order of decreasing 
weight identifies the main wavevectors that characterize the 
state, which can then be checked against the list of optimal 
Luttinger-Tisza wavevectors to help determine whether the 
spin configuration found from Iterative Minimization ought 
to be a ground state. 



3. Common-Origin plot 

As in Ref. 1 , we use the Common-Origin plot to easily see 
which directions the spins point in spin space. To construct 
a Common-Origin plot using a spin configuration determined 
by Iterative Minimization, we simply plot all the spin vectors 
with their tails at the origin. This diagnostic lets us easily 
see which directions the spins point in, but it doesn't give any 
information about how the spins are arranged on the lattice. 

One can also make Common-Origin plots using only a sub- 
set of the entire spin configuration. For example, we might 
make a Common-Origin plot using only spins in one of the 
FCC sublattices. To give another example, in cases where 
the lattice sites break up into sublattices according to site en- 
ergy, we can make a Common-Origin plot using only spins 
with a certain site energy. Making Common-Origin plots of 
subsets of spins like these is often very helpful for finding a 
symmetrical basis for spin space to rotate the spin configura- 
tion into. It is much easier to think of possible idealized forms 
for a ground state when the numerical spin configuration ob- 
tained from Iterative Minimization has been rotated into the 
most symmetrical basis for spin space. 



D. Least Squares Fits and Variational Optimization 

Many of the ground states we find in Iterative Minimization 
simulations are quite complicated, owing to the fact that they 
contain small contributions from many wavevectors besides 
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the optimal Luttinger-Tisza wavevectors. In cases like these, 
it is desirable to obtain an approximate parameterization of the 
ground state. This approximate parameterization can then be 
used to obtain a more complete understanding of the state. To 
find these approximate parameterizations we perform a least 
squares fit of the numerical spin configuration in each F.C.C. 
sublattice to sines and cosines of the optimal Luttinger-Tisza 
wavevectors. If the spin configuration is rotated into a sym- 
metrical basis in spin space, then this method usually yields 
approximate parameterizations which reveal the main struc- 
ture of the ground state. 

When the least squares fit is simple enough (i.e. it only in- 
volves one or two wavevectors), we can idealize the approxi- 
mate parameterization into a guess for the closed-form param- 
eterization of the ground state. Usually this guess involves a 
number of unknown parameters (e.g. wavevectors or cone- 
angles) whose values we compute by analytically solving for 
the values of these parameters which minimize the energy of 
the idealized state. This procedure is known as "variational 
optimization" Q]. 

When we present parameterizations of spin configurations, 
we give a formula S a (r), which represents the spin configu- 
ration as a function of position in the a th F.C.C. sublattice. 
We usually give the components of S Q (r) in a Cartesian ba- 
sis {ei,e 2 ,e 3 } for spin space. Some states, however, take 
on their simplest form when written in terms of a basis for 
spin space which rotates as one moves from site to site in real 
space. It is important to remember that the basis vectors for 
spin space have no special orientation relative to the crystal 
lattice, since equation ( |2.1| > is invariant under a uniform rota- 
tion of all the spins. 



E. Projection of Stacked or Columnar States onto Lower 
Dimensional Lattices 

In this section we discuss a method for analyzing stacked 
and columnar states, in which the spin configuration is inde- 
pendent of one or two of the spatial coordinates, respectively, 
by projecting the lattice and interactions down onto an equiv- 
alent one or two-dimensional lattice. 

A stacked state is one in which there is a certain direction 
in real space, say q, such that within every plane normal to q, 
all the spins have the same direction. Then we can project 
the three-dimensional lattice down onto an equivalent one- 
dimensional lattice (as elaborated in Ref. Q] Section V B), 
in which each spin represents an entire plane of spins from 
the original three-dimensional lattice. The three-dimensional 
spin interactions are correspondingly projected to the one- 
dimensional chain: the interaction %j between Sj and another 
spin Sj in the chain, is the sum of all interactions between 
some particular spin that projects to Sj and any spin in the 
plane that projects to Sj. 

A columnar state is a state in which there is a certain direc- 
tion n, such that in any column of sites parallel to n, all the 
spins point in the same direction. In this case, we can project 
the three-dimensional lattice and its interactions down onto an 
equivalent two-dimensional lattice, in which each spin repre- 



sents an entire line of spins in the original three-dimensional 
lattice. Similarly to the stacked case, the interaction between 
two spins Sj and Sj in this two-dimensional lattice is equal to 
the sum of all interactions in three dimensions between one 
spin in the preimage of Sj and all spins in the preimage of Sj . 

If the result of an Iterative Minimization simulation is a 
stacked or columnar state, then we can make some progress 
towards understanding that state by projecting the lattice and 
interactions down onto a lower dimensional lattice. This lower 
dimensional problem is typically easier to analyze, and in 
cases where the lower dimensional lattice is a Bravais lattice, 
we can immediately use the Luttinger-Tisza method to rigor- 
ously justify the ground state. If a stacked or columnar state is 
the ground state of the three-dimensional lattice, then it must 
also be the ground state of the lower dimensional lattice (but 
the converse is not true). 

The mapping of stacked states was introduced in Ref. 0~| 
for the study of conic spiral states, which break up into a fam- 
ily of parallel planes such that the spin direction is uniform 
within each plane. The state can be optimal in three dimen- 
sions only if its projection is optimal on the one-dimensional 
lattice, but the converse is not true. (That is, a state might 
be optimal on the one-dimensional lattice, but in 3D another 
state that is not stacked in that way might have an even lower 
energy.) 

The one-dimensional lattice can have a non-coplanar 
ground state only if it is non-Bravais lattice, i.e. only if the 
mapped sites are inequivalent by translation. Thus, it appears 
that the parallel layers of spins in the three-dimensional lat- 
tice must be unequal. This happens in the Octahedral lattice 
of Ref. [ 1 1 for (100) layers; in the Pyrochlore lattice the (111) 
stacking consists of triangular-lattice linking layers, alternat- 
ing with kagome-lattice layers that have three times as many 
spins. 

In the case of the Pyrochlore lattice, stacked states in which 
q is a (100) direction project onto a one-dimensional Bra- 
vais lattice; columnar states in which n is a (100) direction 
project onto the two-dimensional non-Bravais lattice known 
as the Checkerboard lattice ll20l . States which are stacked 
along a (111) direction project down onto a one-dimensional 
non-Bravais lattice with a basis of two kinds of site, which 
are projected respectively from a triangular lattice layer of 
spins and a Kagome lattice layer of spins. The ground states 
of ( |2.1| i on this two-site chain lattice, with various couplings, 
were studied extensively in Ref.Q] 



III. OVERVIEW OF RESULTS 

In this section we present a brief overview of our results. 
Most of our work was dedicated to characterizing the ground 
states of ( |2.1| i for all possible combinations of first and second 
nearest neighbor interactions, culminating in the J\- J 2 phase 
diagram. We also briefly explored a few ground states that 
include interactions beyond second nearest neighbors. 
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FIG. 2: The J1-J2 Phase Diagram for the Pyrochlore Lattice. It 
shows the regions of the J1-J2 phase diagram in which the Ferro- 
magnetic (F.M.), Kawamura, q = 0, {gOO} Planar Spiral, Double- 
Twist (D.T.), Multiply-Modulated Commensurate Spiral (M.M.C.S.) 
and Cuboctahedral Stack (C.S.) states are found. Solid lines repre- 
sent first order phase transitions and dashed lines represent second 
order phase transitions. The dashed-dotted line along the negative J\ 
axis indicates that there are a vast number of degenerate ground states 
along that line and that the Kawamura and q = states are members 
of that family of degenerate states. We call this a "degenerate" phase 
transition. 
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A. The Ji- J 2 Phase Diagram 

The J1-J2 phase diagram is shown in figure[2] The ground 
states that we find in this region of parameter space are sum- 
marized in table [I] The two largest regions in the phase dia- 
gram are the Ferromagnetic (F.M.) region and the {q00} Pla- 
nar Spiral region. In the region Ji > 0, J2 > the ground 



state is ferromagnetic (section IV A 1. This is true even in the 



case of a pure ferromagnetic J\ interaction or a pure ferro- 
magnetic J2 interaction. The reason for this is that the J\ and 
J2 interactions act between spins in different F.C.C. sublat- 
tices, so they couple together all of the spins in the lattice. 
The Ferromagnetic state extends into the J\ < 0, J2 > and 
J 1 > 0, J2 < regions of the J1-J2 phase diagram, although 
it is quickly destabilized by the antiferromagnetic interactions 
found in those regions. 



The {^00} Planar Spiral (section IV E 1 is a stacked state in 
which the spins spiral in a plane in spin space as one walks 
in a particular (100) direction in the lattice. This state exists 
in both the J\ < 0, < and Ji > 0, J2 < regions of 
the J\- J2 phase diagram. In the special case where J\ = 
and J 2 < 0, the {g00} Planar Spiral is a 120° spiral in a 
(100) direction in the lattice. This special case corresponds to 
the negative J 2 axis in figure [2] which lies within the region 
where the {g00} Planar Spiral is the ground state. 



Ground State 


Ji 


J2 


NCP? (y/n) 


Qlt. 


Ferromagnetic (F.M.) 


+,- 


+,- 


n 


(0,0,0) 


Ji < 


— 





n 


Any 


{g00} Planar Spiral 


+,- 


- 


n 


{g00} 


q = 






n 


(0,0,0) 


Cuboctahedral Stack (C.S.) 


+ 


— 


y 


aiii 
1 2 2 2 / 


Kawamura States 




+ 


y 


{ffo} 


Double-Twist (D.T.) 


+ 




y 


{ffo} 


Multiply-Modulated 
Commensurate Spiral (M.M.C.S.) 


+ 




y 


13 1 11 
1 4 2 4/ 



TABLE I: Summary of the ground states found in the J1-J2 phase 
diagram. It tells the region of the phase diagram where each state is 
found, whether or not the state is non-coplanar (NCP), and also gives 
the optimal Luttinger-Tisza wavevector Ql.t. for each state. 



The next largest region in figure [2] is the region in the 
Ji < 0, J2 > quadrant where the ground states are the 
three Kawamura States (section[Vl|; the Sextuplet-q state, the 
Quadruplet-q state of type 1 and the Quadruplet-q state of type 
2. The Fourier Transforms of these three states consist primar- 
ily of {ffo} wavevectors. At J\ w — 1.09Ji the Kawamura 
States transition to the ferromagnetic state. 

In the upper part of the J\ < 0, J2 < quadrant of fig- 
ure [2] the ground state is the q = state (section IV Ci, in 
which all spins in the same F.C.C. sublattice are parallel to 
each other, so that the spins point in only four directions in 
spin space. These four spin directions also add to zero be- 
cause of the antiferromagnetic J\ interaction . At J 2 = 0.5 Ji 
this state transitions to the {y00} Planar Spiral state. 

Finally, there are three more ground states in the J\ > 0, 
J2 < quadrant of figure [2] The first of these is the Cuboc- 
tahedral Stack (C.S.) state (section|V]i, a non-coplanar state in 
which the spins point towards the eight corners and 12 edge- 
midpoints of a cube in spin space. As J2 is made more nega- 
tive there is a transition to the Multiply-Modulated Commen- 



surate Spiral (M.M.C.S.) state (section Villi, a non-coplanar 
state in which the spins perform two different kinds of spirals 
which are controlled by two distinct wavevectors. Lastly there 
is the Double-Twist (D.T.) state (section [VIlJl, a non-coplanar 
state in which the spins perform two different kinds of spirals 
but those two kinds of spirals are controlled by the same kind 
of wavevector. 

We have arranged our discussion of these states roughly in 
the order of least complicated to most complicated. In the 



section dealing with the simpler states (section IV 1 we have 
grouped together states which are closely related. So, for ex- 
ample, the discussion of the q = state comes directly af- 
ter the discussion of the ground states of the pure J x < 
antiferromagnet, of which it is a special case. The Cubocta- 
hedral Stack state is the most complicated state for which we 
still have an exact parameterization and energy per site, so our 
section on this state (section [V} can be said to bridge the gap 
between the less complicated and more complicated states. 
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B. States with J3 and J3 



B. Pure Ji < 



In our exploration of states involving the third nearest 
neighbor interactions J3 and J3, we again found the Cubocta- 
hedral Stack state, and we also found a new kind of Alternat- 
ing Conic Spiral state. We found that one can also stabilize 
the Cuboctahedral Stack state with a ferromagnetic J\ inter- 
action, a ferromagnetic J3 interaction and an antiferromag- 
netic J3 interaction. We conjecture that the two regions of pa- 
rameter space in which the Cuboctahedral Stack is the ground 
state (i.e. this region and the region in the J\-J<i plane) are 
smoothly connected. 

When the magnitude of the antiferromagnetic J3 interac- 
tion becomes much larger than the magnitudes of the ferro- 
magnetic J 1 and J3 interactions, we find a new kind of Alter- 
nating Conic Spiral state (appendix |BJ. This state is different 
from the alternating conic spirals discovered in Ref. [J in the 
sense that the wavevector which controls the spiraling behav- 
ior and the wavevector which controls the alternating behavior 
are not parallel to each other. 



IV. SIMPLE STATES 

We start our tour of the phase diagram with the most ele- 
mentary cases, in that there is only one coupling, and/or there 
are ground states with the periodicity of the unit cell. All these 
states are exact Luttinger-Tisza states whose exact energies 
are known, and in none of these cases is a non-coplanar state 
forced. 



The ground state of \2.\ \ in the case of a pure antiferromag- 
netic first nearest neighbor interaction J\ < is well known 
to be any state in which the sum of the four spins on each tetra- 
hedron is equal to zero. One can prove this using the Cluster 
method [7|, in which one finds a way to rewrite the Hamil- 
tonian as a sum of disjoint terms: if a state can be exhibited 
in which each term is separately minimized, this must be a 
ground state (and all other ground states must have the same 
property) 

Let the tetrahedra in the lattice be indexed by the Greek 
letter 77 and let the spin in F.C.C. sublattice a on the if h tetra- 
hedron be s r; Q . Then the pure- Ji Hamiltonian takes the form 



n 



constant , 



(43) 



where Li v = J2 a =o s v,°" l e - me sum °f ^ e f our s P ms on the 
rf h tetrahedron. 

Eq. \A3\ is a sum over disjoint terms, each of which is min- 
imized by 







(4.4) 



Furthermore it is easy to find an example configuration that 
satisfies condition \AA\ simultaneously on all tetrahedra, so 
we can rigorously conclude this condition is true for all 77 in 
all ground states. In any such a state, the site energy is £ \ = J\ 
for every spin, so the energy per site is also 



£ = Ji 



(4.5) 



A. Ferromagnetic State 

The simplest state we find (and the simplest possible ground 
state) is the ferromagnetic state, in which every spin points in 
the same direction. This state is obviously the ground state in 
the region where J x and J 2 are both greater than zero, but we 
also find it in the region J\ < 0, J 2 > when J 2 > — 1.09 Ji 

and in the region Ji > 0, J 2 < when J 2 > f _ | + 1?) Ji- 
The energy per site of the ferromagnetic state is 



£ = -3 Ji - 6 J 2 



(4.1) 



If the interactions are Ji, J 3 and J3 (section V D 1, the energy 
per site of the ferromagnetic state is 



-3../i — 3J3 — 3(/o 



(4.2) 



We use these expressions for the energy per site of the fer- 
romagnetic state to predict where this state transitions to com- 
peting states in the regions J\ > 0, J 2 < and J\ < 0, 
J 2 > 0, and also in regions of the phase diagram where J3 
and J3 are non-zero. 



These states have an extensive degeneracy, meaning the num- 
ber of parameters needed to specify the spin configuration is 
proportional to the number of spins in the system. The q = 
states of Sec lIVCl and the Kawamura states of Sec. I VII still 
satisfy the constraint ( |4.4| i so they are specific subsets of the 
ground states for the pure J\ < case. 



C. q = State(s) (Ji < and J 2 < 0) 

In the region J\ < 0, | J\ < J 2 < we find a state in 
which all spins in the same F.C.C. sublattice are parallel and 
the spins from the four F.C.C. sublattices all sum to zero. To 
be more precise, we have S a = n a , where the n a are all 
constant unit vectors which satisfy 



(4.6) 



We call this state a q = state because the wavevector 
characterizing the spin configuration in each F.C.C. sublattice 
is the zero wavevector. The site energy for each spin in this 
state is £j = Ji + 2 J 2 so this state has an energy per site equal 
to 



S = Ji + 2J 2 



(4.7) 
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For sufficiently small J2, at least, the ground state must be 
a subset of the pure- Ji ground state manifold of Section [TV B| 
selected by J 2 as a degenerate perturbation. This state is most 
quickly understood by using the equivalence of a small anti- 
ferromagnetic (ferromagnetic) J 2 interaction with a small fer- 
romagnetic (antiferromagnetic) J3 interaction, within that de- 
generate manifold (see Ref. [8] for a proof of this fact). The 
reason for this equivalence is the following. The J2 and J3 in- 
teractions of spin % can be gathered into contributions from six 
second-neighbor tetrahedra, each of which includes a differ- 
ent J 1 neighbor of spin i. But in view of ( |4.4[ ) the sum of the 
two J2 neighbors, the one J3 neighbor, and the one J\ neigh- 
bor in each of those tetrahedra is zero, hence the sum of the 
J\, J2, and J3 energies is also zero - but the total J\ energy 
is a fixed constant within the degenerate manifold satisfying 

<@3). 

So, we can trade this problem for that of a small J 3 > 0. 
Since the J3 bonds connect sites of same F.C.C. sublattice, the 
J3 term is optimized when all spins in the same F.C.C. sublat- 
tice are parallel (so q = by definition). This is compatible 
with \\A\ . the necessary and sufficient condition is Eq. ( |4.6| l, 
as was claimed. 

Thus the q = states form a continuous two-dimensional 
manifold of degenerate states, inequivalent by rotational sym- 
metry, and parametrized by two angles, 9 and <fr. We can 
define 9 € [0, ir] to be the angle between n and ri! - this 
must also be the angle between n 2 and 113, since |no + ni | = 
|ri2 + 113 1 ; then we can take <f> £ [0, tt] to be the angle between 
the plane of (no, ni) and that of (112, 113). 

Because the q = family satisfies the constraint < |4.6) , this 
state is also a ground state of j2.l\ in the case of a pure first 
nearest neighbor interaction Ji < 0. Evidently, the set of 
q = states is just a subset of the ground states for the pure 
Ji < interaction, and these states are selected out by the 
antiferromagnetic J 2 interaction. 

Figure [3] shows the four eigenvalues A„ (q) of the matrix 
J(q) plotted along the T^X^W-^L-^T^K^X 
route through the Brillouin zone of the F.C.C. lattice for the 
parameter values Ji = —1, J2 = —0.25. We see that the 
largest eigenvalue has its maximum at the T point q = 0, so 
the optimal L.T wavevector is Ql.t. = 0, not only for small 
J2 but in a considerable interval. Since we constructed nor- 
malized spin states using just these modes, they are rigorously 
the ground states, and Iterative Minimization simulations con- 
firm this. 



D. Pure J 2 < 

Unlike the case of a pure J3 interaction, in which the 
Hamiltonian breaks down into uncoupled sublattices, the pure 
J2 interaction acts between sublattices, so that every spin is 
coupled to every other one. In the case of a pure antiferro- 
magnetic second nearest neighbor interaction, we find that the 
ground state is a 120° planar spiral stacked in a (100) direction 
in the lattice. We can get an understanding of this state by as- 
suming (based on evidence from Iterative Minimization sim- 
ulations) that the state is stacked along a (100) direction and 



FIG. 3: The four eigenvalues A„(q) of the matrix J(q) plotted along 
the T^-X^W^L^T^K^X route through the the 
Brillouin zone of the F.C.C. lattice for Ji = — 1, J2 = —0.25, 
where the ground state is the q = state. Recall that V = (0, 0, 0), 
X = 2tt(1,0,0), W = 2tt(1, |,0), K = 2tt(§, §,0) and L = 
27r(|, |, |). The largest eigenvalue of J(q) takes on its maximum 
value at the F point q = 0, confirming that the optimal L.T. wavevec- 
tor in this region is Ql.t. = 0, as expected. 




r x w l r k x 



then projecting the interactions down onto an equivalent one- 
dimensional lattice. We can then apply the Luttinger-Tisza 
method [3] |6) to this equivalent one-dimensional lattice. We 
use the method outlined in subsection lll El to determine the ef- 
fective interactions in the equivalent one-dimensional lattice. 

Let us take the stacking to be along the z -direction. The 
(001) planes of sites with the same spin direction are sepa- 
rated by a distance of j. A spin has four of its second nearest 
neighbors in each adjacent plane (separated by a distance 4 
in the z direction) and two more second nearest neighbors in 
planes that are two steps away (i.e. separated by h in the z 
direction). Thus, the equivalent one-dimensional lattice has 
both first and second neighbor interactions, ji = 4J 2 and 
= 2J2 = hji< In addition, the one-dimensional lattice 
is a Bravais lattice, because every site in the same constant-z 
plane has the same number of second nearest neighbors in the 
planes above and below it. 

The optimal wavevector takes the form q = (0, 0, q) in the 
three-dimensional lattice. In the one-dimensional lattice the 
state will be characterized by a wavenumber qid- We have 
Q = Qid/% since in the three-dimensional lattice the mini- 
mum separation in the z direction between two spins in the 
same F.C.C. sublattice with the same x and y coordinates is 
|, whereas the minimum separation between two spins in the 
one-dimensional Bravais sublattice is ~. Since the minimum 
separation in the z direction between two spins in the same 
sublattice is doubled when going from the one-dimensional 
lattice to the three-dimensional lattice, the wavevector is 
halved. 

Applying the Luttinger-Tisza method to the one- 
dimensional chain gives an expression for J(qid), the 
Fourier transform of the couplings in the one-dimensional 



9 



chain: 



J(qi d ) = 4 J 2 ( 2 cos 



(qid\ 



(4.8) 



Maximizing this we find that the optimal (one-dimensional) 
wavenumber is given by 



fqid\ 



(4.9) 



so we have qid = which means that spins separated by a 
distance of 4 in the z-direction (first nearest neighbors in the 
one-dimensional lattice) are rotated 120° from each other. So 
the ground state is a 120° spiral in the z-direction. This state 
has the form 



Sid(z) 



'8tt 

f(» | — Z | C| 



'8tt 

sin I — z I c 2 



(4.10) 



in terms of the interactions in the three-dimensional lattice. 
The parameterization of this state in the three-dimensional 
lattice has the same form as equations ( |4.1 1 a| > and ( |4.1 lb} , 
but with ^ replaced by the wavenumber q determined from 
( |4.14| > and the relation q = qid/2- The site energy for each 
spin in this state is 



4J 2 



3J 2 



(4.15) 



so £ = £i. In the special case of J\ = 0, we recover the 120° 
spiral found in the previous section. 

Finally, we note that in order for equation ( |4.14| i to have a 
solution for qid, we must have — 1 < — Jl ^j 2 J2 < 1. There- 
fore, this state can only exist in the region J\ < when 
Ji < \J\, and it can only exist in the region J\ > when 
J 2 < -iJi. 



in the one-dimensional lattice, up to an arbitrary (constant) 
phase angle. In the three-dimensional lattice we have q = ^p, 
and we can parameterize this same state as 



S Q =o,3(r) 
S Q= i, 2 (r) 











= cos 






— sin 










= cos 


— z 




)* 




V 3 


+ 3, 














+ 


sin 













4tt 



ei4.Ha) 



e 2 (4.11b) 



in the separate F.C.C. sublattices. 

The site energy of each spin is given by = 3 J 2 so this 
state has an energy per site 



V. CUBOCTAHEDRAL STACK 

In the region Ji > 0, J 2 < we find a state that we 
have named the Cuboctahedral Stack because the spins in 
each Kagome lattice layer perpendicular to a certain (111) di- 
rection in real space are arranged in the Cuboctahedral state 
found on the Kagome lattice in Ref. |25] and on the Octahe- 
dral lattice in Ref. [T| In each of those states the spins point 
towards the 12 vertices of a cuboctahedron or, equivalently, 
the 12 edge-midpoints of a cube. If we choose a basis in 
spin space in which the sides of this cube are perpendicular 
to (100) directions, the spins point in the twelve (110) direc- 
tions in spin space. 



£ = 3J 2 



(4.12) 



The 120° spiral is just a special case of the more general 
{q00} planar spiral state, discussed in the next subsection 
( |IV E| i, which involves both J\ and J 2 interactions. 



E. {§00} Planar Spiral (J 2 < and J ± > or J a < 0) 

In the region J\ < 0, J 2 < \j\ and also in the region 
J\ > 0, J 2 < — 0.68 Ji, we find a planar spiral state which 
is stacked in a (100) direction, which (as in Sec. |IV D i we 
again take along z. This generalized spiral can also be un- 
derstood by mapping the interactions down onto an equivalent 
one-dimensional lattice, now with couplings j\ = 2Ji + 4 J 2 
and j 2 = 2J 2 . Applying the Luttinger-Tisza method to this 
one-dimensional lattice gives 



J(qu) = 2ji cos + 2j 2 cos v } , 



( qu\ 



(4.13) 



so the optimal one-dimensional wavenumber is given by 

(qid\ ji Ji + 2J 2 

ms = = 



V 4 J 



4j 2 



4J 2 



(4.14) 



A. Structure of the Cuboctahedral Stack 

In the Cuboctahedral Stack state there is a distinguished di- 
rection in real space, which is one of the four (111) direc- 
tions (we choose [111] for the parameterizations below). In 
the Kagome lattice layers stacked perpendicular to this direc- 
tion, the spins point towards the 12 edge-midpoints of a cube. 
In the triangular lattice layers between the Kagome layers, the 
spins point towards the eight corners of that same cube. 

This state is built out of three of the four {^|^} type 
wavevectors and the fourth unused wavevector points in the 
stacking direction of the Kagome and triangular lattice layers; 
it can be written as 

1 3 

S = — = cos(q fe • r) e fe (5.1a) 

v 3 k=l 

1 3 

S a =i,2,3 = —p- ^2 cos(q fc • r) e fc , (5.1b) 



where the three wavevectors used in this parameterization 

are qi = 27r(-i, \, §), q 2 = 2tt{\,-\,\) and q 3 



n i\ 1 _1\ 
Z7 H 2 ' 2' 2>- 
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In this parameterization, the triangular layers consist of the 
spins in F.C.C. sublattice 0, whose configuration is given by 
( |5.1a[ ), and the Kagome layers consist of the spins in F.C.C. 
sublattices 1, 2 and 3, whose configurations are given by 
( j5.1b| l. The site energies for spins in the four F.C.C. sublat- 
tices are 

£[ 0) = -V6Ji (5.2a) 
£^ 2 ' 3) = ~Ji (l + ^ , (5.2b) 

where the superscript indexes the F.C.C. sublattices. The site 
energy of a spin in a triangular layer is lower than the site 
energy of a spin in a Kagome layer. 

B. Energy per spin in the cuboctahedral stack 

The energy per site in this state is then 

Evidently, the J2 contribution to the energy of this state com- 
pletely cancels out, so that the total energy depends only on 

In fact, the lower bound on the energy of this state (from 
equation ( |2. 10| >) is — 2J\, so even though the Cuboctahedral 
Stack is constructed solely from the optimal L.T wavevec- 
tors, it does not achieve the lower bound on the energy. This 
is related to the fact that only three out of a total of four avail- 
able symmetry related {555} modes are used to construct 
this state. 

Figure [4] shows the four eigenvalues A„ (q) of the matrix 
J(q) plotted along the T^X^W^L^T^K 
— > X route through the Brillouin zone of the F.C.C. lattice 
for the parameter values J\ = 1, J 2 = —0.2. We see 
that the largest eigenvalue has its maximum at the L point 
q = 27r(|, i, i), confirming that the optimal L.T. wavevector 
is Q l.t. — 27r(|, |, |), as expected from Iterative Minimiza- 
tion simulations. 



1. Why energy is independent of J2 

To see why the total energy ( |5.3[ ) is independent of J2, we 
need to think more about the structure of this state. Of the 
four spins in every tetrahedron, three of them lie in a Kagome 
lattice layer and the fourth lies in a triangular lattice layer. 
Recall that every spin in a triangular lattice layer points in a 
(111) direction and every spin in a Kagome lattice layer points 
in a (110) direction. To approximately satisfy the ferromag- 
netic first nearest neighbor interaction, the three spins on a 
tetrahedron in a Kagome lattice layer will point in three (110) 
directions which surround the same corner of a cube in spin 
space. For example, these spins might point in the [110], [101] 



FIG. 4: The four eigenvalues A„(q) of the matrix J(q) plotted 
along the T — > X — > W — > L ^ F — > K — > X route through 
the Brillouin zone of the F.C.C. lattice for Ji = 1, J 2 = —0.2, 
where the Cuboctahedral Stack (C.S) state is the ground state. Re- 
call that T = (0,0,0), X = 2tt(1,0,0), W = 2ir(l,|,0), 
K = 27r(|, f , 0) and L = 27r(§, \, \). The largest eigenvalue of 
J(q) takes on its maximum value at the L point, q = 27t(k, 5, h), 
the wavevector which characterizes the Cubocatahedral Stack state. 
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and [Oil] directions. Then the dot product between any two 
of these three spins will be | . The fourth spin on this tetrahe- 
dron, the one which lies in a triangular lattice layer, will then 
point in the [111] direction, the average of the directions the 
other three spins point in. The dot product of this fourth spin 
with any of the other three will then be ^ ps 0.82, which also 
favors the ferromagnetic first nearest neighbor interaction. 

Because of this structure, the dot product between any spin 
in a triangular lattice layer (i.e. in F.C.C. sublattice 0) and any 
of its second nearest neighbors is always zero. The reason for 
this is that any second nearest neighbor of a spin in a triangular 
lattice layer will be in a Kagome lattice layer, so it will point 
in a (110) direction. Moreover, the (110) direction it points in 
will be one which surrounds a corner of the cube in spin space 
which is adjacent to the corner that the spin in the triangular 
lattice layer points towards, but this (110) direction will not 
point towards the midpoint of the edge that connects these 
two corners of the cube in spin space. It follows from this that 
the dot products of a spin in a triangular layer with its second 
nearest neighbors are all zero. 

The dot products between any spin in a Kagome lattice 
layer and its 12 second nearest neighbors are not all zero but 
they do add to zero. Of these 12 dot products, four of them 
are (the dot products with spins in the adjacent triangular 
layers), four of them are — I (the dot products with spins in 
the same Kagome layer) and four of them are \ (the dot prod- 
ucts with spins in the neighboring Kagome layers which are 
beyond the adjacent triangular layers). 
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C. Cuboctahedral stack in the ( Ji , J2) phase diagram 

As explained in the last subsection, when we fix J\ and vary 
J2, the energy of the Cuboctahedral Stack state is unchang- 
ing, but the energy of any competing states (whose energy 
generally does depend on J 2 ) may cross lower than that of the 
Cuboctahedral Stack, at which point there will be a phase tran- 
sition. In particular, comparing Eq. ( |5.3| l to Eq. ( |4.1| >, we find 
that the transition from the ferromagnetic state to the Cuboc- 
tahedral Stack should occur at 




-0.171Ji 



(5.4) 



and this result is confirmed by simulations. Simulations also 
show that the Cuboctahedral Stack transitions to the Multiply- 



Modulated Commensurate Spiral state (section Villi near 
J 2 = -0.40Ji. 



1. Transition between the Cuboctahedral Stack and ferromagnetic 

states 

Because there are no variable parameters in the parameteri- 
zation of the Cuboctahedral Stack state, we might expect that 
the transition between this state and the ferromagnetic state 
would be of first order. The results of Iterative Minimization 
simulations show that this is not the case. On the ferromag- 
netic side of the phase boundary we find a "mixed" state near 
the transition point whose energy is lower than the energies of 
both the ferromagnetic and Cuboctahedral Stack states. The 
Fourier Transform of this mixed state shows peaks at q = 
as well as at the { 5 1 5} wavevectors that make up the Cuboc- 
tahedral Stack. In spin space the spins are arranged around a 
cone so that this state does have a net magnetic moment which 
points along the axis of that cone. 

If one looks in detail at the spin directions in the Common- 
Origin plot, one sees that the spins unfold from the conic axis 
in four sets of three spins and four sets of two spins, with a 
set of two spins between every set of three spins. In addition, 
this configuration is symmetric under a rotation of 90° about 
the conic axis. It appears that as one approaches the phase 
boundary from the ferromagnetic side the four sets of three 
spins become the 12 spins which point in (110) directions in 
the Cuboctahedral Stack state, and the four sets of two spins 
become the eight spins which point in (111) directions in the 
Cuboctahedral Stack state. 

This evidence suggests that the phase transition between 
the ferromagnetic and Cuboctahedral Stack states is continu- 
ous (of second order) and that the spin configuration near the 
phase boundary (on the ferromagnetic side) is a configuration 
which smoothly interpolates between the ferromagnetic and 
Cuboctahedral Stack states. 



D. Cuboctahedral Stack in the Ji, J3, J' A Phase Diagram 

The Cuboctahedral Stack state also exists in the J\, J 3 , J3 
phase diagram for J\, J3 > and J' 3 < 0. This state is 



again made up of three of the four { 5 ^ 5 } type wavevectors, 
with the fourth unused wavevector pointing in the stacking 
direction. The parameterization is exactly the same as that 
presented in equations ( |5.1a| > and ( |5.1b[ ). The energy per site 
for the Cuboctahedral Stack state with these three interactions 



is 




Ji - Ja 



4 



(5.5) 



This state is expected to transition to the ferromagnetic 
state, with energy per site (j4.2|i, when 



Ji 



8J3 



I6J3 



9-2V6 



(5.6) 



From our experience with the phase transition between the 
ferromagnetic and Cuboctahedral Stack states in the J\-J% 
phase diagram, we expect that this phase transition will also 
be of second order. 

Previously, we had understood the Cuboctahedral Stack as 
being the result of competition between a strong ferromag- 
netic nearest neighbor interaction J\ and a weak antiferro- 
magnetic second nearest neighbor interaction J 2 . In this case 
we have no second nearest neighbor interaction, but we do 
have a ferromagnetic J3 interaction and an antiferromagnetic 
J3 interaction. In the Cuboctahedral Stack configuration, a 
spin in a triangular lattice layer has a dot product of | with 
its six J3 neighbors and a dot product of — | with its six J 3 
neighbors. A spin in a Kagome lattice layer has a dot prod- 
uct of 1 with two of its six J 3 neighbors (those in adjacent 
Kagome lattice layers) and a dot product of with the other 
four (those in the same Kagome lattice layer) and it also has 
a dot product of — 1 with two of its J 3 neighbors (those in the 
same Kagome lattice layer) and with the other four (those 
in adjacent Kagome lattice layers). In this way the Cuboc- 
tahedral Stack state serves as a compromise in which spins 
have positive (or zero) dot products with their J 3 neighbors 
and negative (or zero) dot products with their J3 neighbors. 
So the dot product between each pair of spins is either zero or 
has the sign that goes with the interaction between that pair of 
spins (e.g. dot products between J3 neighbors are either zero 
or positive, which favors the ferromagnetic J3 interaction). 



VI. KAWAMURA STATES 

In the region Ji < 0, < J 2 < 1.09 |Ji| of the phase 
diagram, we encounter a remarkable phase first identified in 
finite-temperature simulations of Refs.[36]and[37] and further 
studied in Refs. [8] and [9] This is actually a family of phases 
dominated by a certain kind of {§§0} ordering mode, which 
according to the Luttinger-Tisza analysis is optimal in this pa- 
rameter range. In fact, the exact Luttinger-Tisza wavevector 
starts out at {§§0} for J 2 < |Ji|, but it drifts slightly in- 
wards as J 2 is increased, eventually settling at an incommen- 
surate wavevector near 27r(0.73, 0.73, 0) as J 2 is increased to 
J 2 w \Ji\. So the {||0} wavevectors are only nearly opti- 
mal in this region of parameter space, although they are very 
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close to the true L.T. wavevector. By taking linear combina- 
tions with multiple symmetry-related modes it is possible to 
construct three particularly symmetric noncoplanar states |9| 
that are called the "Sextuplet-q" state, and the "Quadruplet-q 
states of types 1 and 2". We have called this whole family the 
"Kawamura States" after the group which has done the most 
to elucidate their structure. 

The Sextuplet-q state is composed from equal amounts 
of all six {||0} wavevectors. The Quadruplet-q state of 
type 1 is composed from equal amounts of four of the 
{f |0} wavevectors which lie on the points of a cross in re- 



ciprocal space. [For example, these might be 27r(j,4,0), 
27r(f,-f,0), 27r(§,0,§) and 2tt(|, 0, -f ).] Finally, the 
Quadruplet-q state of type 2 is composed of unequal amounts 
of those same four wavevectors. [For example, the largest 
contribution might be from wavevectors 27r(|,|,0) and 
2tt(|, — 1,0), with a slightly smaller contribution from 
27r(|,0, |) and 27r(|, 0, -§)•] Furthermore, each of these 
three states very nearly satisfies the tetrahedron constraint 
(|4.4), so these states are a subset of the pure J\ < anti- 



ferromagnet of section IV B 



Ref. |9] used mean-field theory to predict the stability of 
these phases as a function of temperature, concentrating on 
the particular point J2/IJ1I = 0.2, and found that the cu- 
bic Sextuplet-q phase is stable at higher temperatures, while 
the quadruplet-q kind of state was stable at lower tempera- 
tures. 091 In our T = study, the three states were nearly 
degenerate in this region of parameter space. The L.T. phase 
diagram indicates that the Kawamura states should transition 
to the ferromagnetic state at J2 « 1.09 |Ji|. Iterative Mini- 
mization simulations confirm that a phase transition does take 
place at this location and the simulations also seem to show 
that the transition is of second-order. 

Our Iterative Minimization simulations in this region of the 
J1-J2 phase diagram find these three states, as well as states 
with slightly higher energies which are also mainly composed 
of {|f 0} wavevectors. Therefore we suspect that the energy 
landscape in this region of parameter space consists of three 
local minima representing the three Kawamura states, and that 
these minima sit in a valley with a gently sloping floor. This 
would explain why our simulations sometimes don't find these 
three local minima and instead settle into states with slightly 
higher energies. 

Even though the three Kawamura states are nearly degen- 
erate in this region of the phase diagram, our Iterative Min- 
imization simulations show that the true ground state is the 
Sextuplet-q state for J2 <C \Ji\ at least until J2 = 0.1|Ji|. 
Somewhere between J2 = 0. 1 1 J\ and J2 = 0.15| J\\ there is 
a phase transition, and we find that the true ground state for 
0.15|Ji| < J 2 < 1.09| Ji\ is the Quadruplet-q state of type 
2. It appears that the Quadruplet-q state of type 1 is never the 
true ground state in this region, even though it is always nearly 
degenerate with the other two Kawamura states. 

We have focused our study of the structure of the states in 
this region almost exclusively on the Kawamura Sextuplet- 
q state. There are two reasons for this. Firstly, it has the 
most symmetry in real space and in spin space out of the three 
Kawamura states. More importantly, it has a slightly lower 



Sublattice 


Si 


# of spins 


1 


-1.1053 


192 


2 


-1.1841 


192 


3 


-1.2103 


192 


4 


-1.2672 


64 


5 


-1.2753 


192 


6 


-1.2902 


192 



TABLE II: Site energies and number of spins in each sublattice of the 
Kawamura Sextuplet-q State with Ji = — 1, J2 = 0.1, on a lattice 
of size 4x4x4. 



energy than the other two Kawamura States when J 2 <C | J\\, 
suggesting that it is the state selected out of the highly degen- 
erate manifold of pure-Ji ground states by an infinitesimal 
ferromagnetic J2 interaction. 



A. Structure of Kawamura Sextuplet-q State 

The Kawamura Sextuplet-q state is a complicated non- 
coplanar state which is mostly composed of equal amounts 
of the six { 1 10} wavevectors. Because of the complexity of 
this state, in this section we only present the results of analy- 
sis on a numerical spin configuration generated by an Iterative 
Minimization simulation with the parameter values J\ = — 1, 
J 2 = 0.1 on a lattice of size 4x4x4 with periodic boundary 
conditions. 

In this state, it turns out that the site energy £i of the spins 
takes on only six different values: thus we can organize the 
spins by site energy into six sublattices in real space, each of 
which has the periodicity of a 2 x 2 x 2 supercell (and is not 
itself a Bravais lattice). In fact, the sites in each sublattice are 
equivalent by symmetries of the ground state. The site energy 



the 



and number of spins in each sublattice is given in table 
site positions for each sublattice are given in Appendix 

Sublattices 3 and 4 are the sublattices of highest symmetry 
in this configuration. They are each made up of tetrahedra of 
spins. The tetrahedra in sublattice 4 sit on the sites of a body- 
centered cubic (B.C.C.) lattice with a unit cell of size 2x2x2. 
The tetrahedra of spins in sublattice 3 combine with those in 
sublattice 4 to create a simple cubic lattice of tetrahedra with 
a unit cell of size lxlxl. The spins in sublattice 3 point 
towards the 12 edge-midpoints of a cube in spin space and 
the spins in sublattice 4 point towards the eight corners of that 
same cube. 

The nearest neighbors of the spins in sublattices 3 and 4 are 
distributed throughout the six sublattices. Because sublattices 
3 and 4 are made up of tetrahedra, three of the six nearest 
neighbors of each spin in these two sublattices lie in the same 
sublattice as that spin. The other three nearest neighbors of 
the spins in sublattice 4 (there are 192 = 3x64 of them) are 
contained in sublattice 5 (these are the only spins in sublattice 
5). The other three nearest neighbors of each of the spins in 
sublattice 3 (there are 576 = 3x192 of them) are distributed 
throughout sublattices 1, 2 and 6. 
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Q 


A a 


B a 


C a 


D a 


E a 


F a 


1 


—a 


C2 


-Cl 


-C 2 


Cl 


-c 2 


2 


-Cl 


C2 


C2 


Cl 


-C 2 


Cl 


3 


C2 


— Cl 


— Cl 


— ca 


-c 2 


Cl 


4 


C2 


-Cl 


C2 


Cl 


Cl 


-C 2 



TABLE III: Coefficients for the approximate parameterization ^6.2\ 
of the Kawamura Sextuplet-q state. The greek letter a indexes the 
four F.C.C. sublattices and the coefficients are all given in terms of 
the two constants ci ~ 0.73 and C2 ~ 0.27. 



Next we present the results of a least-squares fit of this nu- 
merical spin configuration to functions of the form Cj sin(q 3 • 
r + <pj) where q,, is one of the wavevectors of type {§§0} 
(an optimal L.T. wavevector) and the Cj and <pj are the unde- 
termined fitting parameters. We use the condensed notation 



$™ = 2w [ -x + ~y 



/3 3 . 
$ m , = 2ir[-x--y), 



(6.1a) 
(6.1b) 



and a similar notation for phase functions involving x and z 
or involving y and z, to simplify the presentation of this state. 
Idealizing the results of the least square fit gives an approxi- 
mate parameterization of this state in each F.C.C. sublattice of 
the form 



B a sin($ y2 )| 
D a sm($ xz + j 

4 



ei 
e 2 



+ {c q sin (<3 

+ {E a sin ($ xy + |) + F a sin ($ xy + |) } e 3 (6.2) 
where the values of the coefficients A a , B a , C a , D a , E a and 



F a are given for each sublattice in table III 

We can see from \6.2\ that if one fixes the x and y coordi- 
nates and moves in the z-direction, then the spin configuration 
looks roughly like a superposition of two distorted conic spi- 
rals about the z-axis. A similar statement holds if one fixes 
the x and z coordinates and moves in the y-direction or fixes 
the y and z coordinates and moves in the x-direction. 

The energy per site for this state is £ = —1.2164. 



B. Anharmonic selection as a consequence of normalization 



Although the parameterization ( |6.2[ ) looks relatively clean, 
it is not normalized on all of the lattice sites. This is because 
the Kawamura Sextuplet-q state is not constructed solely from 
{§§0} wavevectors. A closer look at the Fourier Transform 
of the numerical spin configuration found in Iterative Mini- 
mization simulations shows that it also contains small contri- 
butions from wavevectors like 27r(|, |, |), 2ir(\, \, — f ) and 
2-7r(i, |, — i), even though these are not optimal Luttinger- 
Tisza wavevectors. The first two of these wavevectors are 



equivalent by a symmetry of the F.C.C. lattice (reflection 
through the x-y plane), so they give the same eigenvalues for 
the matrix J(q). The third kind of wavevector is not equiva- 
lent to the first two. 

It turns out that wavevectors of these three types are ac- 
tually just linear combinations of three of the optimal L.T. 
wavevectors. For example, we can write 27r(|, — |,Q) + 
27r(f,-|,0) + 27r(|,0, 1) = 2tt(§, -§,§). But we may 
add to this wavevector a reciprocal lattice vector of the form 
27r(— 2, 2, 0) to map it back to a wavevector in the first 
Brillouin zone. We see then that this linear combination 
of optimal L.T. wavevectors is equivalent to the wavevector 
27t(t, hi §)■ It turns out that the small contributions from 
these three kinds of wavevectors can be understood by looking 
at what happens to our approximate parameterization when 
we normalize it. 

To normalize the parameterization (|6.2), we would compute 



Si = 



(6.3) 



for every spin. We can rewrite the magnitude of the spin Sj as 
Si | 2 — 1 ) . Taylor expanding the recipro- 



= 4/1 



cal of this quantity about the point |s;| =1 (where the state 
is normalized), we find 



1 



1 



(6.4) 



So the first order correction to the non-unit vector spin Sj in- 
cludes a term | | 2 s^. It is from this cubic term that we get 
corrections to the state that involve wavevectors which are 
a linear combination of three of the optimal Luttinger-Tisza 
wavevectors (if we wrote out as an exponential Fourier se- 
ries, we would see the addition of wavevectors happening in 
the exponents). 



VII. A NEW VARIETY OF DOUBLE-TWIST STATE 

In the region Ji > 0, —.68 < J2 < —.43, we find a state 
which is chiefly composed of a dominant set of two of the 
six { 1 10} wavevectors and a sub-dominant set of two of the 
other wavevectors of this type, and these four wavevectors all 
lie on the points of a cross in reciprocal space. For example, 
this state might be made from the wavevectors 27r(|, |,0), 
2tt(§ , - 1, 0), 27r(|, 0, |) and 2tt(| , 0, -|) with the Fourier 
amplitudes of the first two wavevectors being approximately 
1.40 times the Fourier amplitudes of the second two wavevec- 
tors. Luttinger-Tisza analysis confirms that wavevectors of the 
type { 1 10} are the optimal wavevectors for this state. 

Although this state is mostly constructed from {ffO} 
wavevectors, it also contains small contributions from { | j0} 
wavevectors. These extra wavevectors aid in normalizing 
this state, since one cannot construct a normalized state from 
{f f 0} wavevectors alone. 

This state shares three characteristics in common with the 
Double -Twist state found on the Octahedral lattice in Ref.Q] 



(1) There is a distinguished direction m in real space and 
also a distinguished direction S3 in spin space. This 
state can be written using a basis for spin space that uses 
the fixed basis vector S3 and two rotating basis vectors 
A(r • m) and B(r • m) which lie in the plane perpen- 
dicular to S3, so that A(r ■ m) x B(r • m) = S3. As 
one moves in the m direction in real space the rotating 
basis vectors A(r • m) and B(r • m) spiral about the 
fixed vector S3 in spin space. 

(2) There is a second direction n in real space orthogonal 
to the m direction, with the property that the spins also 
spiral as one moves in the n direction, although not 
about the S3 axis in spin space. Since the rotating basis 
vectors stay put as one moves in the n direction, this 
second kind of spiral is completely independent from 
the first kind. 

(3) Both distinct kinds of spiraling behavior are controlled 
by the same type of wavevector. 

To summarize, this state gets its name from the fact that the 
spins are tracing out two different kinds of spirals in two di- 
rections which are perpendicular to each other, but these two 
spirals are controlled by the same type of wavevector. In the 
Double-Twist state that we find on the Pyrochlore lattice the 
distinguished direction m in real space is a (100) direction, 
which we take to be the z-direction in our discussion of this 
state. The second spiraling direction n is then the x- or y- 
direction. Both kinds of spiral are controlled by a {f00} 
wavevector. 

We now present an approximate parameterization of this 
state obtained by performing a least squares fit of a numeri- 
cal spin configuration found in an Iterative Minimization sim- 
ulation onto sines and cosines of q • r, where qis a {ffo} 
wavevector (an optimal L.T wavevector). This simulation was 
performed with the parameter values J x = 1 and J 2 = —0.6 
on a lattice of size 4x4x4 using periodic boundary condi- 
tions. 

The two dominant {||0} wavevectors making up this 
state are 27r(|,|,0) and 27r(|,-|,0) and the two sub- 
dominant wavevectors making up this state are 27r(|,0, |) 
and 2?r(f ,0, — |). As in ((64} we use the condensed notation 

<S>x V = 2w(^x+^ (7.1a) 

$ xy = 2tt^x-^ (7.1b) 

$ x = 27r^a:J (7.1c) 

with similar notations for arguments involving y and z. We 
express this state in terms of the two rotating basis vectors 

A» = cos($ 3 )ei+sin($ z )e2 (7.2a) 
BO) = - sin ($ 2 )ei+ cos ($ 2 )S 2 , (7.2b) 
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A' 




C' a 
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di 


d 2 


-dx 


d 2 


d 3 


d± 


5.0° 


3.8° 


2 


di 


d 2 


-di 


d 2 


di 


d 3 


31.0° 


4.5° 


3 


d 2 


di 


-da 


di 


di 


d 3 


17.7° 


-7.4° 


4 


d 2 


di 


-d 2 


dj 


d 3 


di 


18.4° 


17.3° 



TABLE IV: Coefficients and phase angles for the approximate pa- 
rameterization \73\ of the Double-Twist state. The greek letter a 
indexes the four F.C.C. sublattices and the coefficients are all given 
in terms of the four constants di ~ 0.23, d 2 ~ 0.76, ^3 « 0.26, 
di w 0.75. Analysis of multiple examples of this state seems to in- 
dicate that di 7^ d 3 and d 2 7^ di, even though their values appear to 
be very similar. 



which have been chosen so that A(z) x B(z) = S3. The 
approximate parameterization of this state in each F.C.C. sub- 
lattice has the form 

+ {C' a cos ($„ + 6 a ) + D' a cos ($„ + ^ a )} B(z) 

+ [E' a sin (® xy + |) + F' a sin ($ xy + |) } S 3 (7.3) 

where the values of the coefficients A' a , B' a , C' a , D' a , E' a and 
F' a and the phase angles 9 a and t/j a for each sublattice are 
given in table |IV| 

We see from ( |7.3"| that if one holds z constant and moves 
in the y-direction, the spins spiral (but not about the S3 axis, 
since the S3 component of the spins changes as one moves in 
the y-direction), and if one holds y constant and moves in the 
z-direction the spins spiral about the S3 direction (since the 
basis vectors A(z) and B(z) change as one moves in the z- 
direction). In addition, these two spirals are controlled by a 
{ § 00} wavevector. So in this example the m direction is the 
z-direction and the n direction is the y-direction. 

This state has an energy per site £ = —2.2780. 



VIII. MULTIPLY-MODULATED COMMENSURATE 
SPIRALS 

In the small sliver of parameter space Ji = 1, —0.43 < 
J2 < —0.40 we find a non-coplanar state with period 4 
constructed solely from wavevectors of the types (4 \\}, 
{ f|0} and { 4 4 0} (recall that these three types of wavevec- 
tor are not equivalent by any symmetry of the F.C.C. lattice). 
Our prototype of this state uses the four wavevectors q 1+ = 
27r(f,i,i), qi _ = 27r(f, !,-§), q 2 = 2tt(| , -f , 0), and 
q3 = 27r( j, —j, 0), but one can build this state from any 
set of four wavevectors that is equivalent to these by sym- 
metry. The {f f 5} wavevectors are the optimal Luttinger- 
Tisza wavevectors in this region of parameter space. In fact, 
the largest eigenvalues of the matrix J(q) for these wavevec- 
tors at the parameter values J\ = 1, J 2 = —0.42 are 
A maa; (qi±) = A L . T . = 2.0195, A TOaz (q 2 ) = 2.0121 and 
Amaz (q3) = 0.9337. The Luttinger-Tisza wavevectors qi± 
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FIG. 5: Largest eigenvalue X max (q) of the matrix J(q) vs. J2/J1 
for the optimal L.T. wavevector of each of the five ground states 
found in the region Ji > 0, J2 < 0. The name of the corre- 
sponding ground state is shown in parentheses next to the label for 
the type of wavevector. The five ground states are the Ferromag- 
netic (F.M.), Cuboctahedral Stack (C.S.), Multiply-Modulated Com- 
mensurate Spiral (M.M.C.S.), Double-Twist (D.T.) and {q00} Planar 
Spiral states (with wavenumber q for the Planar Spiral given by equa- 
tion < |4. 14| >). The inset shows the small region —0.43 < J2/J1 < 
—0.40 where the largest eigenvalue for the dominant wavevector in 
the Multiply-Modulated Commensurate Spiral state briefly crosses 
above the largest eigenvalue for the wavevectors characterizing the 
Cuboctahedral Stack and Double-Twist states, confirming that there 
is a distinct ground state in this small region of the J\-J% phase dia- 
gram. 
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are the dominant wavevectors in this state, but the state also 
contains sizable contributions from modes containing the sub- 
optimal wavevectors q 2 and q 3 (see equations ( |8.5| l and ( |8.6[ )). 
So this state is another example of a non-coplanar state which 
is constructed using an optimal L.T. wavevector and a set of 
wavevectors which are suboptimal. 

Figure [5] shows the largest eigenvalue X max (q) of the ma- 
trix J(q) for the optimal Luttinger-Tisza wavevector of each 
of the five states found in the region J\ > 0, J2 < as a func- 
tion of J2/J1. It also shows a zoomed in view of the region 
—0.43 < J2/J1 < —0.40 where the largest eigenvalue for 
wavevectors of the type {f 25}, the main type of wavevector 
in the Multiply-Modulated Commensurate Spiral state, briefly 
crosses above the largest eigenvalue of J(q) for the wavevec- 
tors used in the Cuboctahedral Stack and Double-Twist states. 
The fact that the L.T. phase diagram shows a sliver in which 
{III} wavevec tors are optimal is strong evidence that the 
corresponding spin state found in Iterative Minimization sim- 
ulations is not merely an artifact of relaxation. 

We call this state a Multiply-Modulated Commensurate 
Spiral (M.M.C.S.) because the spins in this state are trac- 
ing out two different spirals in two orthogonal directions and 
these two spirals are controlled by two different kinds of 
wavevectors. 



In this state there is a distinguished direction m in real 
space and a distinguished direction S3 in spin space. The dis- 
tinguished direction in real space is parallel to q 1+ + qi_, so 
m points in the [310] direction in our example state. The spin 
configuration can be written using a basis in spin space that 
uses the one fixed basis vector S3 and two rotating basis vec- 
tors A(£) and B(£), where the coordinate £ = r • m. These 
two basis vectors lie in the plane perpendicular to (33 and sat- 
isfy A(£) x B(£) = S3. As one moves in the fn direction 
in real space these two basis vectors rotate counter-clockwise 
about e 3 in spin space. To write this state in the most con- 
cise way we also use a second set of two rotating basis vectors 
A'(£) and B'(£) lying in the plane perpendicular to e 3 and 
also satisfying A'(£) x B'(£) = S3, but this second set ro- 
tates clockwise about S3 as one moves in the fn direction in 
real space. 

There is also a second distinguished direction n in real 
space which is parallel to qi + — qi-,so n points in the [001] 
direction in our example state. Because qi + and qi_ differ 
only in the sign of one component, we have m • n = 0. In this 
state the spins also spiral about about the e 3 direction in spin 
space when one moves in the n direction in real space, even 
though the four basis vectors A(£), B(£), A'(£) and B'(£) 
stay put as one moves in this direction in real space. Finally, 
the wavevectors that control the spiraling in the m and n di- 
rections are not equivalent by symmetry, so the spiraling be- 
havior in this state is different from the spiraling behavior of 
the Double-Twist state of section [VTI| (the two kinds of spirals 
in the Double-Twist state are controlled by the same type of 
wavevector). 

Using a diagnostic which counts the number of different 
spin directions present in a certain state (two spins Sj and Sj 
are considered to point in different directions if s, ■ Sj < .99), 
we find that the spins in this state point in 42 different direc- 
tions. The spins in two of the F.C.C. sublattices point in the 
same 10 directions and the spins in the other two F.C.C. sub- 
lattices point in 32 directions which are distinct from the first 
10. 

Here we present an approximate parameterization of this 
state which was constructed by idealizing the result of a 
least square fit of a numerical spin configuration to sines and 
cosines of q • r, where qis a {§ { ffO}, or a { ~ -0} 
wavevector. This numerical spin configuration was generated 
by an Iterative Minimization simulation with the parameter 
values Ji = 1 and J2 = —0.43 on a lattice of size 4x4x4 
using periodic boundary conditions. 

The approximate parameterization uses the four wavevec- 
tors qi+, qi-, q 2 , and q3 defined in the first paragraph of this 
section. We again use the condensed notation 
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as a shorthand for these position dependent phase angles. The 
most concise presentation of this state uses two sets of two 
rotating basis vectors, all lying in the plane in spin space per- 
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pendicular to (33. We can write these four basis vectors as 

A(x,y) = cos $ xy ei + sin $ xy e 2 (8.4a) 

B(x,y) = - sin $ XJ ,ei + cos (8.4b) 

A'(x,y) = cos $ xy ei - sin $ xy e 2 (8.4c) 

B'(x,y) = sin^ei +cos$ xy e 2 (8.4d) 

where we have arranged it so that A(x,y) x B(x,y) = 
A'(x, y) x B'(x, y) = A x B = e 3 . A' (a;, y) is the reflec- 
tion of A(x, y) over the e! axis and B'(x, y) is the reflection 
of B(x, y) over the e 2 axis. This parameterization takes quite 
different forms in the different F.C.C. sublattices. In sublat- 
tices and 3, we have 

S Q =o,3 = h |sin(vrz)A'(x, y) =F cos(ttz)B'(x, y)| 

+ l 2 |± sin(7rz)A(x, y) — cos(7rz)B(x, y)| 
+ {Z 3 sin6 X2/ + Z 4 sin *.,.,,} e 3 , (8.5) 

where the top (bottom) sign in the plus/minus and minus/plus 
signs corresponds to the zeroth (third) F.C.C. sublattice and 
the coefficients are h sa 0.868, l 2 « 0.054, l 3 sa 0.699 and 
I4 sa 0.076. These are the two sublattices in which the spins 
point in 32 different directions, as discussed above. In sublat- 
tices 1 and 2, the spin configuration has the same form, 

S a =i,2 = 2Z 3 cos(7rz)B'(x,2/) 

+ {l 5 sm@ xy +/ 6 sin* xy }e 3 , (8.6) 

where Z5 sw 0.228, Iq sa 0.011 and Z 3 is the same coefficient 
that appeared in ( |8.5| l. These are the two sublattices in which 
the spins point in the same 10 directions. 

In this parameterization m points in the [310] direction 
and n points in the [001] direction. Furthermore, the spi- 
raling in the m direction is controlled by the wavevector 
27r(|, |,0) and the spiraling in the n direction is controlled 
by the wavevector 27r(0, 0, |), and these two wavevectors are 
clearly not of the same type. 

The energy per site for this state is £ = —2.0131. 



IX. CONCLUSION 

In conclusion, we have surveyed the full range of the phase 
diagram of classical ground states with nearest and second- 
nearest neighbor exchange couplings in the Pyrochlore lattice 
(Figure [2] and Table [Iji. Large swathes of the phase diagram 
(roughly those in which \J 2 \ > | Ji|) are dominated by the 
ferromagnet or a relatively simple coplanar spiral phase, but 
with smaller \J 2 \ values (for either sign of Ji), a zoo of com- 
plex non-coplanar phases were found. 

In a similar study of the Octahedral lattice, which formed 
the model for this one, three categories of non-coplanar state 
were identified, with only one or two examples of each on 
the Octahedral lattice, and our results on the Pyrochlore fit 
into the same categories. One generic category is an ideal 



Luttinger-Tisza state, in which three degenerate ordering vec- 
tors are matched with three spin directions to form a highly 
symmetric, already normalized state with no admixture of 
other modes needed. On the Octahedral lattice these were 
the two cuboctahedral states; our sole example on the Py- 
rochlore is also cuboctahedral (Sec. [V}. The second category 
was multiple-wavevector commensurate states, dominated by 
a star of symmetry-related wavevectors but admixing others. 
On the Octahedral lattice, the only example was the "double- 
twist"; on the Pyrochlore, we have found three representa- 
tives: the Kawamura states (Sec.lVIb, a kind of Double-Twist 



state (Sec. VII 1, and a "Multiply-Modulated Commensurate 
Spiral" (M.M.C.S.) state (Sec. |VIII[>. The third category was 



"conic spirals" which were stackings of layers with a uniform 
spin direction in each layer, which are built from two unrelated 
but nearly degenerate spin modes, one of them being generi- 
cally incommensurate. Being layered, these states can be opti- 
mized after projecting all interactions onto a one-dimensional 
"chain lattice" HI, as we outlined in Sec. II E A corollary is 



that if a non-coplanar state is forced, the layers must be in- 
equivalent, which is indeed the case for the (111) layering in 
which we did find a conic spiral state (Appendix [B]) in the 
Pyrochlore lattice. In both the Octahedral lattice and now in 
the Pyrochlore lattice, it was difficult to stabilize a conic spi- 
ral using the first few neighbor couplings: in the Pyrochlore 
J 3 and J3 were needed. 

When we examine the whole phase diagram, there seems to 
be a rough tendency that the wavevectors used for the (Ji,J 2 ) 
ground state are the same as those used for {—J\, ~J 2 )', for 
example, the Double-Twist state is diametrically opposite the 
Kawamura state in Figure [2] and both use optimal wavevectors 
of type { § |0}; the q = antiferromagnet is opposite the fer- 
romagnet which is (of course) also a q = state. There is a 
plausible reason for this trend. Notice first that the matrix of 
Jij 's (in real space) has no diagonal terms, so its trace must be 
zero. The same must be true for its Fourier transform, so the 
sum of the four eigenvalues of J(q) is zero at every wavevec- 
tor. The optimum wavevector Ql.t. occurs at the point in 
the zone where one branch of the eigenvalue spectrum has its 
greatest positive excursion. Necessarily, the sum of the other 
three eigenvalues must have its greatest negative excursion at 
the same wavevector: not uncommonly, that will be the point 
of the single most negative excursion. But if we reverse the 
signs of all couplings, the L.T matrix is the same except for a 
global sign, and the greatest negative excursion becomes the 
greatest positive one, i.e. the optimal mode. 
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Appendix A: Fourier Transform of the Interactions for the 
Pyrochlore Lattice 

In this appendix we present the matrix J(q) for the Py- 
rochlore lattice including interactions J\, J2, J3 and J3 (i.e. 
including both kinds of third nearest neighbor interaction). 
One can calculate the elements of this matrix using equation 



To concisely express the matrices that go with the two kinds 
of third nearest neighbor interaction, we introduce the con- 
densed notation 

C xy = cos hz±*L \ (A4a) 
C xy = cosf^p^ (A4b) 



M 3 = 



M 3 = 



There are a few interesting things about this matrix. While 
J(q) is generally a hermitian matrix, it is symmetric in our 
case because every site on the Pyrochlore lattice is a center of 
inversion symmetry. Secondly, only J3 and J3 terms appear 
on the diagonals since these are the only interactions that act 
between spins in the same F.C.C. sublattice. The J\ and J2 
interactions only act between sites in different F.C.C. sublat- 



( |2.7| i. The matrix is: 

J(q) = JiMa + J2M2 + J3M3 + J 3 Mjj (Al) 

where 



(A2) 



(A3) 
I 

and a similar notation for cosine terms with arguments using 
q y and q z . 



(A5) 

(A6) 
I 

tices, so these show up only in off-diagonal terms. 

Appendix B: A New Kind of Alternating Conic Spiral State 

Although no conic spiral-type ground states are found in the 
J1-J2 phase diagram of ( |2.1| i, conic spiral ground states can 
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be stabilized on the Pyrochlore lattice for certain sets of inter- 
actions. With the interactions J\ > 0, J3 > and J3 < 0, 
we find a new type of Alternating Conic Spiral state when 
the magnitude of J3 is much larger than the magnitudes of J\ 
and J3. Our state is distinct from the Alternating Conic Spi- 
ral found in Ref. 1 because the wavevector which controls the 
spiraling behavior is not parallel to the wavevector which con- 
trols the alternating behavior. This means that the component 
of the spins perpendicular to the conic axis spiral about that 
axis when you move in one direction in the lattice, and the 
component of the spins parallel to the conic axis alternates in 
sign when you move in a completely different direction in the 
lattice. 

This state is made up of two different wavevectors, qi and 
q 2 . The first wavevector qx controls the spiraling behavior 
and it is likely that qi does not have to be commensurate with 
the lattice (as is the case for the spiraling wavevector in the 
{gOO} planar spiral of section IV El. The second wavevector 
q2 controls the alternating behavior and in order for this state 
to be normalized it must be commensurate with the lattice and 
satisfy cos(q2 • r) = ±1 on all lattice sites r. As mentioned 
above, c\x and q 2 are not parallel. 

We can express this state most concisely using one rotating 
basis vector. If we define the phase $ = qx • r, then we can 
write this basis vector as 



A(r) = cos($)ex + sin(<i>)e2 



(Bl) 



With this rotating basis vector the parameterization of this 
state takes the form 

s = sin a A(r) + cos a cos(q2 • r) e 3 (B2a) 

si = A(r) (B2b) 

S2 = sina A(r) + cosacos(q2 • r) S3 (B2c) 

s 3 = cos(q 2 ■ r) e 3 . (B2d) 



composed of an optimal L.T wavevector, q2, and a second 
type of wavevector, qi, which is only slightly sub-optimal. 
This provides further evidence for the hypothesis advanced in 
Ref.Q]that many non-coplanar states are made from combina- 
tions of optimal Luttinger-Tisza wavevectors and wavevectors 
that are only slightly suboptimal. 



Appendix C: Unit cells of the six sublattices of the Kawamura 
Sextuplet-q State 

These are the locations of lattice sites in the 2x2x2 unit 
cells of the six sublattices for the example of the Kawamura 



Sextuplet-q state example presented above (section VI A 1. 
This spin configuration was generated from an Iterative Min- 
imization simulation with parameter values J\ = — 1, J2 = 
0.1 on a lattice of size 4x4x4. 

1. F.C.C.O: (0,0,0), (1,1,1), (.5,0,1.5), (1.5,1,-5), (.5,1.5,0), 
(1.5,.5,1) 

F.C.C.l: (.5,1.25,.75), (1.5,.25,1.75), (0,.75,.75), 

(1,1.75,1.75), (.5,.75,1.25), (1.5,1.75,-25) 

F.C.C.2: (.25,1,.25), (1.25,0,1.25), (.75,1,1.75), 

(1.75,0,.75), (.25,5,1.75), (1.25,1.5,-75) 

F.C.C.3: (.25,.25,1), (1.25,1.25,0), (.25,1.75,5), 

(1.25,.75,1.5), (.75,1.75,1), (1.75,.75,0) 

2. F.C.C.O: (0,1,0), (1,1,0), (1.5,0,5), (1.5,0,1.5), (.5,5,1), 
(5,1.5,1) 

F.C.C.l: (5,.25,1.75), (.5,1.25,1.75), (0,1.75,.75), 

(1,1.75,.75), (1.5,.75,.25), (1.5,-75,1.25) 

F.C.C.2: (.25,0,.25), (.25,0,1.25), (.75,1,-75), 

(1.75,1,-75), (1.25,5,1.75), (1.25,1.5,1.75) 

F.C.C.3: (1.25,-25,1), (1.25,1.25,1),(.25,.75,.5), 

(.25,75,1.5), (.75,1.75,0), (1.75,1.75,0) 



We can see from ( |B2| i that the spins in FCC. sublattices 
and 2 are performing Alternating Conic Spirals, the spins 
in FCC. sublattice 1 are performing a planar spiral about the 
conic axis, and the spins in FCC sublattice 3 are alternating 
pointing parallel and antiparallel to the axis of the conic spiral. 

An Iterative Minimization simulation at the parameter val- 
ues Ji = 1, J3 = 1 and J3 = —4 on a lattice of size 
4x4x4 with periodic boundary conditions found this new 
kind of Alternating Conic Spiral with qi = 2ir(— |, — |, \), 
q2 = 2n(— i, i, |) and cone angle a s» 0.92 (radians). As J3 
is varied in this region, Fourier Transforms of the numerical 
spin configurations found in Iterative Minimization simula- 
tions still show large peaks near {§§5} type wavevectors. 
For this reason, we suspect that the vector q2 is always a 
{M^} tyP e wavevec t° r (and therefore commensurate with 
the lattice). On the other hand, it appears that qi and the cone 
angle a can vary continuously. 

For this set of parameter values, the optimal Luttinger-Tisza 
wavevector is of the type { 1 1 \ }, so q2 is equal to Ql.t. ■ The 
Luttinger-Tisza eigenvalue is Al t. = 7. We find, however, 
that A maa; (qi) = 6.8629, which is only slightly less than 7. 
So this is another example of a non-coplanar state which is 



3. F.C.C.O: (0,5,5), (0,5,1.5), (0,1.5,5), (1,5,1.5), 
(1,1.5,-5), (1,1.5,1.5) 

F.C.C.l: (0,.25,.25), (0,.25,1.25), (0,1.25,.25), 

(1,.25,1.25), (1,1.25,-25), (1,1.25,1.25) 

F.C.C.2: (75,5,1.25), (75,1.5,.25), (.75,1.5,1.25), 

(1.75,.5,.25), (175,5,1.25), (1.75,1.5,55) 

F.C.C.3: (.75,55,1.5), (75,1.25,5), (.75,1.25,1.5), 

(175,.25,.5), (1.75,55,1.5), (1.75,155,5) 



F.C.C.O 
F.C.C.l 
F.C.C.2 
F.C.C.3 



(0,1.5,1.5), (1,5,5) 
(0,1.25,1.25), (1,55,55) 
(75,5,55), (1.75,1.5,1.25) 
(75,55,5), (1.75,155,1.5) 



F.C.C.O: ((0,1,1), (1,0,0), (.5,0,5), (1.5,1,1.5), (.5,5,0), 
(1.5,1.5,1) 

F.C.C.l: (.5,55,75), (1.5,1.25,1.75), (0,1.75,1.75), 

(1,75,75), (.5,75,55), (1.5,1.75,1.25) 

F.C.C.2: (.25,1,1.25), (1.25,0,55), (75,0,75), 

(1.75,1,1.75), (.25,1.5,1.75), (155,5,75) 

F.C.C.3: (.25,1.25,1), (155,55,0), (.25,1.75,1.5), 

(155,75,5), (75,75,0), (1.75,1.75,1) 
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6. F.C.C.O: (0,0,1), (1,0,1). (-5,1.5), (.5,1,1.5), (1.5,.5,0), 
(1.5,1.5,0) 

F.C.C.l: (1.5,.25,.75), (1.5,1.25,-75), (0,.75,1.75), 

(1,.75,1.75), (.5,1.75,-25), (.5,1.75,1.25) 

F.C.C.2: (1.25,1,.25), (1.25,1,1.25), (.75,0,1.75), 



(1.75,0,1.75), (.25,.5,.75), (.25,1.5,-75) 

F.C.C.3: (.25,.25,0), (.25,1.25,0), (1.25,1.75,-5), 

(1.25,1.75,1.5), (.75,-75,1), (1.75,-75,1) 
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